Canonical form of Hamiltonian matrices 
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On the basis of shell model simulations, it is conjectured that the Lanczos construction at 
fixed quantum numbers defines — within fluctuations and behaviour very near the origin — smooth 
canonical matrices whose forms depend on the rank of the Hamiltonian, dimensionality of the vector 
space, and second and third moments. A framework emerges that amounts to a general Anderson 
model capable of dealing with ground state properties and strength functions. The smooth forms 
imply binomial level densities. A simplified approach to canonical thermodynamics is proposed. 
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Since many physical problems can be thought as ma- 
trix diagonalizations, much attention has been devoted to 
the general properties of matrices, starting with Wigner's 
Gaussian Orthogonal Ensemble (GOE) which has pro- 
vided invaluable insight into fluctuation properties [|l|,|[. 
However, physical systems are usually associated to 
Hamiltonians of rank (r) 1-1-2 (i. e., l-|-2-body) while 
for GOE the rank is the dimensionality of the space. 
Two Body Random Ensemble (TBRE) or Embedded 
GOE (EGOE) are the suggested denominations for lim- 
ited ranks. We shall prefer to speak of Hamiltonian ma- 
trices, either random or not. Exact results and even nu- 
merical ones are hampered by the difflculty of enforcing 
the rank conditions, which arc hard to simulate through 
special band or sparse matrices, and it would be desirable 
to have a canonical form that ensures rank 1-1-2. It is our 
aim to show that the Lanczos tridiagonal construction 
can provide it. (The impatient reader may proceed to 
the discussion after Fig. |l|.) 

The algorithm is standard |^ : 

A starting vector (the pivot) |0) is chosen. The vec- 
tor \ao) = HlO) has necessarily the form |ao) = Hoo\0) + 
|1') >, with (0|1') = 0. Calculate (0|ao) = i^oo- Normal- 
ize |1'), to find Hoi = (I'll')^/^. Iterate until state \k) 
has been found. The vector \ak) — 'H\k) has necessarily 
the form 



\ak) = Hkk-i\k - 1) + Hkk\k) + \{k + 1)'), 



(1) 



Calculate {k\ak) = Hkk- Extract and normalize |(fc-f 1)'), 
to find Hfefc+i = ((fc + l)'|(fc + i)r/2. 

The choice of pivot is arbitrary and it can be adapted 
to special problems. One of the most interesting is the 
calculation of strength functions: Given a transition op- 
erator T, act with it on a target state \t) to define a pivot 
|0) = T\t) that exhausts the sum rule for T. The tridiag- 
onal matrix elements are then linear combinations of the 
moments of the strength distribution (0|7i''^'|0) (expand 
|0) in eigenstates to check the statement). Therefore the 
eigensolutions of the / x / matrix define an approximate 
strength function Si{E) = Y.^=l,I ^(E - Ei){i\T\0)^, 
whose first 2/ — 1 moments are the exact ones (state |/) 



is obtained by orthogonalizing T-i^^^\0) to the preceeding 
states. Hence Hjj = (0|7i^^~^|0) + information from 
previous iterations) The eigenstates act as "doorways" 
whose strength will be split until they become exact so- 
lutions for / large enough. To fix ideas consider in Fig. n^ 
the first 200 Lanczos matrix elements calculated to re- 
produce the Gamow Teller strength (T = ctt) obtained 
through the reaction ''^Ca(p,n)''^Sc. In Ref g 700 itera- 
tions proved necessary to ensure the exact S{E) in the pf 
shell (the dimension of the JT = 13 space is d = 8590). 
The question at the origin of the present work is: do 
we truly need 700 exact eigensolutions to obtain a good 
S{E)? Probably not, and the idea is to calculate ex- 
actly a few matrix elements, and replace the rest by an 
adequate approximation. 
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FIG. 1. Exact matrix elements for ^*Sc (J — 1) 



The very rapid variation of the first matrix elements 
in Fig. n certainly depends on the pivot. It is followed 
by fluctuations around smooth patterns: a constant for 
the diagonals, and a line b for the off diagonals. There- 
fore, a good approximation to the exact problem comes in 
three steps. First determine the smooth patterns, then 
do the few Lanczos iterations necessary to reach them 
and finally replace the rest of the matrix elements by 
fluctuations mounted on the smooth patterns. 



Next we concentrate on the first step. 

The quantities that turn out to be necessary are the 
dimensionality, which we write as d = 2^, the centroid 
Ec = d^^ti-{H), the variance cr^ and skewness coefficient 

71 : 



rf = 2^ =^ iV = In d/ In 2, E^ = d'hriH) 

a^ = d-hiiH - Ecf, 71 = {da^)-^lv{H - E^f . 

An inverse binomial function is also needed: 
X = nib(&, TV), < a; < 1/2, the inverse of 

b = bin(.T) = e^[(l+2^)ln(l+2^) + (l~2:c)ln(l-2a;)]^ 

the Stirling approximation of 



(2) 
(3) 



(4) 



h = bin(x, N) = 2 



-N 



N 
N{l/2~x) 



a binomial shifted and normalized so as to be maximum 
and unity at the origin. 

The "data" from which we draw conclusions are ex- 
act diagonalizations (at fixed angular momentum J) for 
"'^Ca, a typical closed shell nucleus (s), and for the config- 
urations {f7/2P3/2)t{99/2d5/2)t ( [fpgdf for short) that 
exhibit rotational behaviour with backbend [|6[ . 

The matrix elements in Figs. and correspond to 
maximum d for each family. They are totally represen- 
tative of the other members. Only the general (pivot- 
independent) trend is visible. 

To within fluctuations the diagonals are constant or 
logarithmic, and the off-diagonals have an inverse bino- 
mial form. 

The matrix elements that reproduce — perfectly within 
fluctuations — the exact ones are 

F,, = -7i(ln^ + l), H,,+^ = ^'^/Nmh{-,,N). (5) 
a d 

The subtraction in the In term ensures that Hi i is trace- 
less, a' and 7J are determined so as to yield a and 71 
through 



a^ - {re) - d-^ Y. ^'+1 + ^'-1 + H. 



2 

a ' 



(6) 
Hi)- (7) 



For the symmetric case (71 w 0), the solution is simply 
a' ~ a. For the asymmetric case, Eqs. (O and M) have 
to be solved. All parameters are given in the captions to 
Figs. I and I 

An approximation to the nib function is Gauss^^ ~ 
X = yJ—{2N)~^ In g , the inverse oi g — exp(— 2A^a:^) ob- 
tained by expanding h in Eq. (0) around a; = 0. Therefore 
Hii^i in Eq. (phis approximated by 



which is seen in Figs, g and || to be fairly close to 
nib except for the severe overshooting at the origin that 
will become worse as N increases. To fix ideas, think 
of the quotient of the direct functions: at x = 1/2, 
5(l/2)/6(l/2) = 2^e'^/2 « e°-2^. 
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200 400 600 800 1000 1200 1400 1600 
FIG. 2. Exact matrix elements for **Ca (J = 4), the in- 
verse binomial nib= Ha+i (Eqs. (0)) , and Gauss"^ = H^^j 
(Eq.( 1)). Parameters: d = 1755, A^ = 10.78, a = 7.85 MeV 
7i — —.01. See text 







v/-(a')V21n(z/d), 



(8) 
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FIG. 3. Exact matrix elements (only one every ten are plot- 
ted) for ifpgd)^ (J = 6) compared with log= Ha (arbitrarily 
shifted), and nib= Ha+i from Eqs. (bl), and Gauss^^ = H^^i 
(Eq. (H)) . Parameters: d = 6579, iV = 12.68, a' = 2.14 MeV, 
a = 2.18, 7( = -0.246, 71 = -0.296. See text 



Eqs. (|h are not fits but checked conjectures that 
should lead to the following CAT (Computer Assisted 
Theorem): 

CAT. Hamiltonian matrices with no hidden symme- 
tries have, except near the origin and within fluctuations, 
canonical Lanczos forms that depend only on their di- 
mensionality, second and third moments and rank of the 
interaction. 

No hidden symmetries imply that the matrices cannot 
be split in blocks. Which means that eigenvalues can- 
not be degenerate 0, page 300]. In other terms: the 



spectrum must have level repulsion. If this condition is 
not fulfilled, degeneracies (not detected by the Lanczos 
method), or semi-degeneracies (leading to strong fluctua- 
tions) will make it impossible to define a canonical form. 
Therefore, it is here that the condition of good quantum 
numbers comes in. 

The rank dependence is crucial; e. g., for GOE matri- 
ces the equivalent of Eqs. <^ is Ha = and Ha+i — 
a/1 - i/d {sic). (The cherry on the CAT cake) 

The behaviour of the matrix at the origin depends on 
the pivot. The only general statement is that Ha+i must 
be (O = of order) 0{a) at the origin and raise to reach 
the nib line after 0{N) iterations. If the first diagonal 
is taken to be the lowest unperturbed state, or chosen 
variationally, it will be at an energy 0{N) below the 
log line and raise to reach it after 0{N) iterations. It 
seems a rather general fact that for such "good" pivots, 
diagonals and off-diagonals reach their respective lines 
in « N/2 iterations Rpl. It goes without saying that 
the correct descriptions of low-lying states and strength 
functions hinges on a good description at the origin. The 
contribution of Eqs. (^) is nonetheless essential: trun- 
cated calculations will have smaller a and the nib line 
will be lower. As a consequence, the ground state energy 
will not converge to the correct value and the density 
of states that determine the final aspect of the strength 
functions will be badly missed [g| . 

Eqs. (|5|) — eventually supplemented by the rapid raise 
at the origin — define an integrable problem, leading to a 
spectrum of locally equidistant levels (picket-fence) and 
locally uniform strength distributions. When fluctua- 
tions are switched on, the level spacings can be expected 
to follow a Wigner law, and the strengths a Porter- 
Thomas distribution [gj. To put it in other words: the 
system becomes chaotic and the strength is localized. 

Eqs. (H), duly "dressed", amount to a general but con- 
strained Anderson model, in which constants are replaced 
by log and nib functions, while behaviour at the ori- 
gin and fluctuations (the dress) must be characterized 
in terms of properties of the Hamiltonian. It is likely 
that the large oscillations, 0{N) steps near the origin — 
as in Fig. ^ — depend both on the pivot and on "details" 
of the Hamiltonian. Then, secular trends take over: fluc- 
tuations in the diagonals are about double that of the 
off-diagonals, and both seem to follow a random walk 
in which — at each step — it is more probable to keep the 
same direction than to reverse it. 

Perfect dressing amounts to exact solutions, but model 
dressing is also possible and can be quite instructive. 

For level densities, fluctuations matter less and it is ob- 
vious that the digonalization of the matrices defined by 
Eqs. (P) will yield densities very close to the exact ones. 
It is perhaps not so obvious that CAT has a corollary: 

CAT corollary Level densities at fixed quantum num- 
bers are binomial. 



The checks are fully convincing for all the calculated 
cases , except when the dimensionalities become very 
small {d ~ 100). Fig. is an example. 
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E (MeV) 
FIG. 4. Exact level densities for (fpgd)^, J = 6 compared 
with the binomial ones. Parameters from Fig. H, treated as 
explained in [hll Eq. 5] (do = 1) 



The only significant problem is that the position of the 
ground state — which depends critically on the behaviour 
of the matrix at the origin — is not predicted with suffi- 
cient accuracy for spectroscopic purposes. 

It should be noted that in Fig. |j an Edgeworth cor- 
rected Gaussian Q works as well as the binomial pO[ . 
As explained after Eq. (g), binomials and Gaussians — 
and a fortiori Edgeworth corrected ones — can be close 
for the small N accessible to simulations. But not in 
general. 
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FIG. 5. Level densities for [fpgdf, J = 0, 12. The legend 

indicates the order in which pj appear as dj increases (full 

lines) or decreases (dashes) 



For (fpgd)^, the calculated multimodal total density, 
Pm — d^^ X](2'^+ 1)pji turns out to be as well described 
by a single binomial as the unimodal ones. The gen- 
tle buildup shown in Fig. H helps in understanding what 
happens. For ''^Ca the situation is the same. We have 



two examples of a multimodal pm behaving as unimodal, 
which squares with the systematic success of Bethe's for- 
mula jl^] . Assuming that the signal for an eventual phase 
transitions is a significant deviation from binomial be- 
haviour, the conclusion is that no low energy thermal 
phase transitions occur in nuclei. 

For the individual pj, Bethe's expression O, Eq.(2B- 
62)] amounts to Ec = J{J + 1) and constant ct^, while — 
for ifpgd)^ — in Fig. ||, we find that both Ec and tr^ 
are nearly parabolic, the former with opposite sign to 
Bethe's. For ^*Ca, Ec is nearly constant and cr^ again an 
inverted parabola. No reason for alarm, but the subject 
demands further study. 



partition function for a general problem is enormously 
simplified, as it becomes the sum over x, a set of con- 
served quantum numbers (Eq^ is the energy origin): 




5 10 15 20 25 

FIG. 6. Centroid E^, c^^ -671 -2 and d/1000 for (fpgdf. 
Full (open) squares for J odd (even). 

One thing is clear: we have to learn how to calculate 
low moments at fixed quantum numbers. A brute force 
approach is viable in very large spaces, but a bit of for- 
malism could extend the range to arbitrarily large ones. 
The dimensionalities in Fig. |g provide an example: they 
are the differences of two binomials, obtained by realizing 
that the calculation of the number of states at fixed Jz 
is a single particle problem. 

Once we have the moments, many things can be done. 
One of them is to treat random Hamiltonians and un- 
derstand the canny insistence of J = to be ground 
state Q. Probably the most promising is the calcula- 
tion of thermodynamic quantities. 

Thermodynamically, the behaviour of a single binomial 
is elementary. Write p{n) = A"(jJ, with E = ne. Then 
z(/3) = EnPHe''^"" == (1 + Ae-'3^)^, is the partition 
function and E{f3) = d\nz/d(3 = A^eAe-''V(H- Ae-'^^), 
the thermal energy: The same as obtained from a micro- 
canonical calculation, /3 = d\np/dE, using the Stirhng 
approximation for p. Therefore, continuous and discrete 
binomials have the same thermodynamics; The canonical 



^(/5)=^(l + A,e-'3^==)^^e 



/3£ox 



(9) 



More on this soon |l4 



Let us close by making explicit the striking comple- 
mentarity between the approach in the previous paper, 
involving scalar traces d~^{'H^), and the Lanczos ap- 
proach here, involving expectation values for a pivot at 
fixed quantum numbers, (0|7i^|0). Low values of K de- 
scribe the centroid region, and are easy to find for the 
traces, while they describe the ground state in Lanczos 
and are hard to calculate. The converse is — or at least 
seems to be — true for high K. CAT is more general and 
has far more consequences than the CLT results [O , but 
they are theorems, while CAT is a conjecture. The theo- 
rems will be of help if they can be extended to traces at 
fixed quantum numbers. A task we have to tackle at any 
rate. 
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